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Abstract 

We use mutual information as a measure of the entanglement between 'physical' and 
thermo-field double degrees of freedom in field theories at finite temperature. We com- 
pute this "thermo-mutual information" in simple toy models: a quantum mechanics 
two-site spin chain, a two dimensional massless fermion, and a two dimensional holo- 
graphic system. In holographic systems, the thermo-mutual information is related to 
minimal surfaces connecting the two disconnected boundaries of an eternal black hole. 
We derive a number of salient features of this thermo-mutual information, including 
that it is UV finite, positive definite and bounded from above by the standard mutual 
information for the thermal ensemble. We relate the construction of the reduced den- 
sity matrices used to define the thermo-mutual information to the Schwinger-Keldysh 
formalism, ensuring that all our objects are well defined in Euclidean and Lorentzian 
signature. 
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1 Introduction 

Mutual information (MI) provides a valuable measure of the entanglement of a quantum 
system. Unlike the von Neumann entanglement entropy of spatial regions, mutual informa- 
tion does not generically suffer UV divergences in quantum field theory (QFT). The mutual 
information between two subregions A and B describes the total correlation, both classical 
and quantum, between the regions [1]; in particular, the mutual information provides an 
upper bound on the correlation for any two observables supported in A and B respectively 
[2]. As a result the mutual information may be used to ascertain basic properties of any 
local observable in a quantum system. 

Unfortunately, a practical method of computing MI for arbitrary regions in a generic 
QFT does not exist, making detailed analysis of MI in systems of interest quite difficult. 
Exact analytic computations of MI are rare and often rely heavily on conformal symmetry 
(see e.g. [3, 4, 5] and references therein). Following the conjecture of Ryu and Takayanagi 
[6] it is widely believed that the AdS/CFT correspondence can be used to compute quantum 
informatics quantities in certain strongly coupled QFTs. In particular, it is believed that 
the Ryu- Takayanagi (RT) formula [6] for the holographic entanglement entropy captures 
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the leading iV 2 contribution the entanglement entropy of the dual CFT. Efforts have been 
made towards a concrete derivation of the RT formula (see, e.g., [6, 7, 8, 9, 10]) but this 
formula and it's covariant formulation [11] remain a conjecture. Furthermore, currently there 
does not exist a tractable method for calculating subleading corrections to the RT formula 
(see for instance [9]). However, the robustness of the RT formula has inspired a great 
deal of investigation into the relationship between quantum information and holography 
[12, 13, 14, 15]. 

In this paper we examine the holographic mutual information between regions on dis- 
connected boundaries in the maximally extended eternal Schwarzschild-AdS spacetime. The 
standard holographic interpretation of asymptotically AdS spacetimes with multiple discon- 
nected boundaries is that they provide a bulk realization of an entangled state in a system 
of multiple non- interacting field theories 1 . In the simplest case of a maximally extended 
Schwarzschild-AdS black hole, which has two causally disconnected boundaries, one can in- 
terpret one boundary as the thermo-field double of the field theory on the second boundary 
[17]. Therefore, the mutual information we study holographically should characterize the 
mutual information between operators and their thermo-field double - we call this type of 
mutual information "thermo-mutual information" so as to distinguish it from mutual infor- 
mation for purely 'physical' operators. The thermo-mutual information quantity provides 
information about the level of entanglement of the thermal system. It is our hope that this 
work will lay necessary foundation for a more broad study of holographic mutual information 
between regions on disconnected boundaries of bulk spacetimes. 

As a warm-up to our analysis of holographic systems, we compute the mutual information 
between 'physical' fields and their thermo-field doubles in two simple systems: a quantum 
mechanical two-site spin chain and a two-dimensional massless fermion. We then calculate 
thermo-mutual information for a 2d holographic system dual to asymptotically AdS^ space- 
time in the simple case of a non-rotating BTZ black hole. We find similar behavior in these 
three very different systems. 

The rest of the paper is organized as follows. Section 2 is devoted to introducing thermo- 
mutual information in the setting of quantum field theory without regard to holographic 
duals. We begin by briefly reviewing mutual information and thermal quantum systems in 
§2.1 and §2.2 respectively; we then define thermo-mutual information and describe it's basic 
attributes in §2.3. We provide two example calculations of TMI: the first in the comfortable 
setting of a quantum mechanical spin chain (§2.3), and the second, after describing some 
computational details in §2.5 (some derivations are in Appendix A for brevity), in a 2D CFT 
(§2.6). Section 3 is devoted to studying holographic systems. We first review the conjecture 
of Ryu and Takayanagi in, and give a brief summary of the non rotating BTZ geometry, 
laying out the tools which make our calculation straightforward. The explicit calculation of 
mutual information and thermo-mutual information is in §3.2. We close with a discussion 
and comments on generalizations of our calculations in section 4. 

1 The robustness of this interpretation is debated [16]. 
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2 Thermo- mutual information in field theory 



Unless stated otherwise, throughout this section we consider QFTs on a space A4 = R x A 
with R the continuous time and A a (d — l)-dimensional spatial lattice with spacing a, which 
we may simply take to be square. The QFT has a Hilbert space H and an algebra of local 
observables A(Ai) represented by bounded operators on H. The reason for considering a 
lattice theory is that for such theories the Hilbert space that is explicitly a direct product of 
the Hilbert space at each site, H = ®hH S itei and therefore may be decomposed into spatial 
regions. The support of a local observable may be as small as a single lattice site labelled 
by the position vector x. We expect to be able to address continuum QFTs by taking the 
a — > limit; more carefully, we could say that we define a continuum QFT by this limit. 

2.1 Mutual information 

Mutual information is constructed from reduced density matrices, so we begin by discussing 
these objects. Quite in general a quantum state \I/ may be described by it's density matrix 
p : = p(\l/; A{M)). The lengthy notation denotes that p is properly defined as a functional 
on A(Ai) which defines the state Correlation functions of an observable O G A(Ai) with 
respect to \I/ may be computed by tracing against p: 

{0) v :=ix H \pO\. (2.1) 

However, if we restrict our attention to observables supported in a subset A C T, t with E t 
an equal-time slice on A4, then we may compute correlation functions without knowing p 
in it's entirety; it suffices to have the reduced density matrix pa := p(ty; A(A)). Due to the 
underlying lattice structure of M. the Hilbert space may be factorized in position space, i.e. 
H = Ha x H A where A is the complement of A on S f . Then pA is obtained by tracing over 
H— ' 

P A = tTnJp}. (2.2) 

Correlation functions of any observable Oa £ A(A) may then be obtained from operations 
in Ha alone via 

(0 A ) 9 = tT HA \p A A ]. (2.3) 
The mutual information for two disjoint subsets A, B C T> t is given by 

Ml{A : B) := S{p A ) + S{p B ) - S(paub) > 0, (2.4) 

where S(p) is the entanglement entropy (a.k.a. Von Neumann entropy) associated to a density 
matrix: 

S(p) := -tr[plogp]. (2.5) 

The trace is over the Hilbert space for which p is a trace operator. The final inequality in 
(2.4) is a consequence of the fact that entanglement entropy satisfies strong subadditivity, 
i.e., 

S(A) + S(B)> S(AUB). (2.6) 
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Mutual information quantifies the total correlation, both classical and quantum, between 
the subsets A and B [1]. This is seen, for instance, in the Pinsker inequality [2]: 

m{A :B)>\ ( {0a °^~^\ {0b) * ) 2 , O a e A(A), O b e A(B). (2.7) 

Mutual information may be defined for continuum QFTs by taking the limit a — > 0. In 
this limit the entanglement entropy S(p) is UV divergent, under reasonable conditions [18] 
scaling as S(pa) ~ Area.(dA)/a d ~ 2 + . . .. For d = 2 QFTS, which will be what we study in 
detail, the divergent part of S(p) has the universal scaling behavior 

S( PA ) ~ 7 ciV log a. (2.8) 

Here 7 is a constant which depends on the theory but not the specifics of the region A, c is 
number which characterized the "number of effective degrees of freedom" (i.e. the central 
charge in a CFT), and N is the number of boundary points of A. As a result of this universal 
behavior we may construct a UV-finite symmetric function [19]: 

F(A : B) = S(p A ) + S(p A ) - S(p AuB ) - S(p AnB ). (2.9) 

The mutual information agrees with this function within Mi's domain of definition Af]B = 0. 
Similar statements a believed to hold generically in higher dimensions as well, as all divergent 
terms are expected to be expressible as integrals over OA that only depend on the cutoff 
scale and therefore should cancel - see discussion in [4, 20, 21] - though these statements 
are tempered by the fact that there exist few exact calculations of entanglement entropy. 

2.2 Thermal systems 

The canonical description of thermal states in QFT follows from the semi-classical inter- 
pretation of these states as equilibrium states in a grand canonical ensemble. Consider a 
system governed by a time-independent Hamiltonian H. The thermal state Q with inverse 
temperature /3 is defined by a density matrix p = e~^ H ; normalized correlations functions of 
the state are constructed via: 

• • .)„ := |tr[p0(a;) . . .] = ^tr[e^ H cj>{x) . . .], (2.10) 

where Z := tr[p] so that (l) n = 1. This description of Q provides a simple connection to 
classical statistical mechanics. 

However, it is often preferable to describe Q as a pure rather than mixed state. This 
may be accomplished with the formalism now known as "thermo-field dynamics" (TFD) - 
original works include [22, 23], for useful reviews see [24, 25]. TFD employs an enlarged 
quantum system composed of two copies of the original Hilbert space %. The total Hilbert 
space of this system is given by the tensor product 

U TFD = U l ®U^ (2.11) 
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with Hi-, I-L2 isomorphic to %. The total Hamiltonian is 

H TFD = H x ® 1 2 - U <g> # 2 , (2.12) 
and f2 may be written in the product space as 

* n 

where \E n ) i are eigenstates of Ifj belonging to %i respectively. From (2.12) and (2.13) we 
readily see that H^fd \Q) = 0. For every operator 4>(x) of the original theory there are 
two copies <pi(x), 4>2{x) in the 'purified' theory. Correlators involving only type-1 operators 
reproduce correlators of the original theory, e.g., 

(n\ Mx)Mv) ■■■\ty = (<Kx)<t>(y) . . .) n . (2.14) 

Since the Hamiltonian does not mix operators on % x an d H2, all type-1 operators commute 
with all type-2 operators: 

(n\...[Mx),Mv)]---\n) = o. (2.15) 

The algebra of observables of the TFD system is 

A T fd(M) := Ax(M) ® A 2 (M). (2.16) 

It is not always conventional to regard type-2 operators as "observable." Because correlators 
of type-1 operators reproduce correlators of the original theory, the second copy of the field 
theory is sometimes referred to as just a 'useful fiction' or computational trick. Furthermore, 
correlators of exclusively type-2 operators may be related to correlators of exclusively type-1 
operators via a anti-unitary map. Ultimately, our motivation for placing the type-2 operators 
on equal footing with those of type-1 originates from the dual description in terms of eternal 
black hole spacetimes where the two types of fields live democratically. 



2.3 Thermo- mutual information 

We can now define the thermo-mutual information (TMI). Consider two subsets A, B C 
and the associated sub-algebras of observables A\{A) G Ai(M) and A2(B) G A2{M). Using 
these ingredients we may construct the reduced density matrices 

p M := p^AM)) (2.17) 
p B2 := p(n;A 2 (B)) (2.18) 
PAWB2 := p(n; A^A) U A 2 (B)). (2.19) 

Since Hi may be factorized in position space we may decompose H.\ = Wai ®^3h e ^c; then 
these density matrices are simply 

Pax = tr« Xl tr Wa p, (2.20) 
PB2 = tr Wl tr w _ 2 p, (2.21) 
Paiub2 = tr H tr H p. (2.22) 
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The first two density matrices correspond to reduced density matrices of the type considered 
in §2.1; the third is a hybrid. We define the thermo-mutual information to be: 



TMI(v4i : B 2 ) := S(p A1 ) + S(p B2 ) - S(p AluB2 ). 



(2.23) 



TMI has the same structure as MI but we reserve the notation MI to refer to the case where 
the two subalgebras Ai(A), Aj(B) belong to the same sector i — j. We could equivalently 
define TMI with the structure of the function F(A : B) (2.9) as the last term S(p A mB2) 
vanishes due to the fact that A\(A) D A 2 (B) = 0. 

Thermo-Mutual information is indeed a kind of mutual information. In particular: 

i) TMI is non-negative: Since *4i(^4) and A 2 (B) are independent their associated den- 
sity matrices satisfy strong subadditivity (2.6) and thus 



ii) TMI bounds local correlators between A and B: This is once again a result of the 
Pinsker inequality. For O A \ £ A\{A) and B 2 £ A 2 (B) the Pinsker inequality yields 
the bound 

( 1 2) " 2 I \\O a1 \\\\Ob2\\ ) • ( 5) 

Other salient features of TMI are: 

iii) TMI provides a lower bound for MI When A and B are disjoint the following 
inequality holds: 



This is a simple consequence of the fact that the tensor product A\{A)®A 2 (B) is always 
equally or more disordered than Ai(A)®A\(B). While we do not have an explicit proof 
of this, it holds in all cases we study, and likely follows simply from subadditivity and 
the explicit form of the thermo-field double state. 

iv) TMI is free of UV singularities: Just as with mutual information and the F function, 
in the limit a — > the entanglement entropies in (2.23) have UV divergences. So long 
as the UV divergences are indeed universal and geometric as described in §2.1, these 
divergences cancel within TMI(Ai : B 2 ) for all subsets A and B, even when A fl B ^ 0. 
To see this note the following two facts. First, the "number of effective degrees of 
freedom" in the type-1 and type-2 sectors are equal. Second, because A(Ai) and A(B 2 ) 
are independent the geometric region relevant to the divergent terms in S(paiub2) is 
vol(cM + dB) rather than vo\{d(A U B)) as in the case of ML 

v) TMI is bounded above for f3 > 0: This follows from (iii) and (iv) provided that the 
theory does not suffer divergences associated to large separations (infrared divergences). 
The j3 — > limit is discussed below. 



TMI(Ai : B 2 ) > 0. 



(2.24) 



MI(Ai : Ei) - TMI(Ax : B 2 ) = S(A 1 U B 2 ) - S(A l U B x ) > 0. 



(2.26) 
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vi) TMI is non-zero for some A, B: This is a trivial corollary of (ii), and is necessary 
for Q to be entangled. Assuming < (3 < oo then Q does not factorize into a product 
state on %tfd Assuming A(M) is non-empty, then there exists an observable O such 
that (0 1 2 ) n > 0. 

We also note the behavior of TMI as (3 — > oo and (3 — > 0: 

vii) TMI vanishes as f3 — > oo (temperature T — > 0): In this limit the canonical ensemble 
becomes a pure state of lowest energy; equivalently, in the TFD formalism Q — > \E ) 1 £g> 
\Eq) 2 . Pure states have maximal correlation, so in this limit the mutual information 
tends to it's maximum value. In contrast, the TMI tends to zero in this limit because 
the type-1 and type-2 sectors are no longer entangled. 

viii) The /3 — > (T — > oo) limit: In this limit the canonical ensemble becomes random 
and there are no (connected, normalized) correlations. For disjoint regions both the MI 
and TMI tend to zero (with the MI bounding the TMI from above). When the regions 
overlap we generically find that the TMI diverges as /? — > 0. This does not contradict 
the points above. 

We may summarize these points by stating that TMI provides a robust measure of the 
correlation between type-1 and type-2 fields. Two important roles of the TMI(Ai : B2) are 
that it bounds from below the MI(A : B), and it bounds from above the correlation between 
any pair of observables in Ai(A) and *4. 2 (Z?). 



2.4 Example in quantum mechanics: the two-spin system 

As a first introduction to TMI let us consider an example in the simple setting of quantum 
mechanics. Consider a system of two "sites" each with a two-spin degree of freedom governed 
by the Hamiltonian 

H = S A ■ S B . (2.27) 

The sites are labelled A and B and the spin vector Sa = S\ + S y A + S Z A and likewise for S^- 
More details than one could ever want to know about this system can be found in any good 
quantum mechanics textbook (e.g., Ch. 3 of [26]). Each site has two states which we label 
by the eigenvalues ± of S A , S B , and so the two-site system has four tensor-product states. 
The energy eigenkets of the system are 2 

= Tsf \+) A ® Hb ~ Ha ® \+)b\ e i = -3, 

\a 2 ) = ^(\+) A ® \~)b + \~)a ® \+)b), E 2 = h (2.28) 
l«3) = |+) A ®|+) B , ' E 3 = l, 

\<H) = |-) fl , E 4 = l. 



2 The energy eigenkets correspond to the simultaneous eigenkets of (Sa + Sb) 2 and (S Z A + SI) 
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We first consider the system as a canonical ensemble. The density matrix describing the 
thermal state Q with inverse temperature /3 is 

-PH ^ a En 4 
P=-^ = J2-^\ a n)M, Z = J2e- f3E " = e^ + 3e^. (2.29) 

n=l n=l 

The entanglement entropy of the entire system is just the thermodynamic entropy which is 
easily computed: 

S(P) = ^ {3 log(3 + e^) - log (l - ^) } . (2.30) 

If we trace of one of the sites then we are left with a random ensemble; the reduced density 
matrix is (in the spin basis) is pa = diag{l/2, 1/2}. The entanglement entropy of a single 
site is therefore 

S( PA ) = S(p B ) = \og2. (2.31) 

From (2.30) and (2.31) we obtain the mutual information MI(A : B) = 2 log 2 — S(p). 

Now let us apply the TFD formalism to the system. The state Q becomes an entangled 
pure state 

I 4 

= -pt £ e ~ PEn ® > (2-32) 
v Z 1 

v 71=1 

To compute the TMI we need the density matrix for A\ U B 2 ; it is a straightforward if 
aggravating exercise to trace over the B\ and A 2 spin sites in order to obtain this density 
matrix. This matrix is not diagonal in the A\ <8> B 2 spin basis, but is easily diagonalized. 
The entanglement entropy may then be computed: 



S (Paiub2) = ~^^\ (-3 + e 2 2 log 



(-3 + e 2 ^) 2 
4(3 + e 



+ 3(l + e 2/3 ) 2 log 



cosh 2 (/3) 



4cosh(2/3) -2sinh(2/3) 



(2.33) 



and from this we obtain TMI(Ai : B 2 ) = 2 log 2 — S(paiub2)- 

We plot MI(Al : B x ) and TMI(Aj, : S 2 ) as a function of /3 in Fig. 1. All features of this 
plot are as anticipated: 

i) The TMI bounds the MI from below; this may be analytically verified by examining the 
quantity 

m(A x : £1) - TMI (At : B 2 ) = -S{p) + S(p AluB2 ) > 0, (2.34) 
with S(p) and S(paiub2) as in (2.31) and (2.33). 
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Figure 1: The mutual information and thermo- mutual information in the two-spin system. 
The solid blue line is the MI and the solid red line is the TMI. The solid black line at the 
top of the graph denotes 2 log 2 which is the maximum value of the mutual information. The 
dashed blue and red lines denote the correlation functions \ {(S A1 S Bl ) n ) 2 and \ {(S A1 S B2 ) n ) 2 
respectively. 



ii) Also plotted in Fig. 1 are the right-hand sides of the Pinsker inequalities (2.7) and (2.25) 
for the observables S Z A1 S Z B1 and S A1 S Z B2 . F° r the case at hand these are simply 

1 z z 2 (e 4/3 - l) 2 1 z 2 2(e 2/3 - l) 2 

2 (\ S Ai S Bi)n) = 2( g 4( 3 + 3)2' 2 (^i^^n) = ( e 4/3 _|_ 3)2 • ( 2 - 35 ) 

These are examples of observables bounded above by the MI and TMI respectively. 

iii) In the (3 ->• limit both the MI and TMI tend to (recall point (viii) of §2.3). 

iv) In the — > oo limit the mutual information tends to it's maximal value of 2 log 2 while 
the TMI tends to zero (recall point (vii) of §2.3). 



2.5 Computing TMI in field theory 

The computation of entanglement entropy in field theories is considerably more involved 
than in the simple quantum mechanical example above. The standard prescription known 
as the "replica trick" involves the following steps (see, e.g., [3]). One first constructs a 
path integral representation of the reduced density matrix p A '■= p( x I / ; A(A)). Using this 
representation one may compute the "Renyi moments" tr^ A (p^). This procedure amounts 



10 



1 


[2- -i 




■ 
■ 
■ 


1 m 


1 1 

1 IF 



Figure 2: Schwinger-Keldysh time integration contours in the complex t plane. The Eu- 
clidean contour is the dashed red line; the TFD contour is the solid blue line. These contours 
are further discussed in Appendix A. 



to gluing together path integral representations such that tr nA (p 1 \) maybe be viewed as a 
single path integral over a Riemann surface rather than the original space M.. Assuming the 
expression for tr^ A (p^) is a complex analytic function of n in a connected region R of the 
complex n plane which includes the non-negative integers, one may obtain the entanglement 
entropy S(p A ) via 

d 

S(pa) = -ti nA p A \ogp A = - lim — tr^ A (p^). (2.36) 

n-+l on 

The path integrals involved in this prescription are on most solid footing in Euclidean sig- 
nature. 

In order to utilize the replica trick to construct both the MI and TMI we need a path 
integral representation of the reduced density matrices p A i, Pbi, Pb2, Paiubi, and Pavob2- 
Life is simpler if we work in Euclidean signature Ai — > A^e = S 1 x A where the Euclidean 
time circle has radius 0. The only subtlety that arises is understanding how type-1 and 
type-2 fields are distinguished in Euclidean signature. Recall that the Euclidean and TFD 
formulations of thermal field theory are related by a deformation of the Schwinger-Keldysh 
time contour - see Fig. 2. From this perspective there are two types of operators in Lorentz 
signature because the Lorentzian time contour has two legs. However, there is just one leg to 
the Euclidean time contour and therefore only one set of Euclidean operators. So care must 
be taken to define the Euclidean analogue of a density matrix in which one type of Lorentzian 
operators has been traced over. At the end of the day, the Euclidean path integrals defining 
these density matrices are an intuitive generalization of existing results [3, 4], so we simply 
quote the results below and relegate a careful derivation to Appendix A. 

We use as a coordinate on Me the complex coordinate z with Rez = x G K and 
Im z G (—(3,0). For A a single interval the density matrix p A \ has the Euclidean path 
integral representation 



PAi(tp',^) 



[V<j>]e 



-S E [. 



II 5 (V( 0+ ) - ^(O + )W0((T) - ^(o-) 



(2.37) 



Here denote Euclidean field profiles and we have labelled only the dependence on 

Imz, i.e., 0(O + ) = 0(Rez,Imz = + ie), etc. Away from the delta function insertions the 
Euclidean fields <f)(x) satisfy the KMS boundary conditions <p(z) ~ (f)(z — i/3). Eq. (2.37) 
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is precisely the path integral obtained from the zero-temperature case by compactifying the 
Euclidean time direction [4]. The expression for paiubi follows suit: 



Paivbi{iP',^) = Z 1 J [T>4>]e 



■SeI 



J] sU(o + ) - ^'(0 + ))5U(0-) - ^(o-) 



n 5(V(0+) - ^'(0+))5(V(0-) - V(O-) 



(2.38) 



The path integral for a type-2 density matrix differs from (2.37) in two respects: the delta 
function insertions are placed at Imz = —/3/2, and the ie prescriptions are reversed: 



p A3 W,1>) = Z- 1 [V<j>)e 



Y[s(4>(-f3/2 + 0-)-^(-f3/2 + 0- 



x£A 



5(0(-/3/2 + O + )-^(-/3/2 + O + )) 



(2.39) 



Of course, using a coordinate transformation one may recast (2.39) in the form of (2.37), so 
the Renyi moments of pai and p^2 are equivalent. These differences do matter, however, in 
the hybrid density matrix: 



PAiuB2&',i>) = Z- 1 / {V<p]e 



n 5(0(0+) -^( O + )j^0(O-)-^(O 

xEA 

J] s((j>{-0/2 + 0") - il/{-i0/2 + 0-)) 
res 

sL{-iP/2 + 0+) - i/;(-ip/2 + 0+)) 



(2.40) 



The placement of delta insertions is depicted in Fig. 3. 

Given (2.37)-(2.40) it is straightforward to proceed with the remaining steps of the replica 
trick. The result is simply that we study the partition function on a surface constructed by 
gluing thermal cylinders together, the only difference is that not all cuts are at Imz = 0, 
and the orientation of the cuts at Imz = —(3/2 are reversed. This is a consequence of our ie 
prescription, and is very similar to calculations of entanglement negativity, where one acts 
with a partial transpose on a reduced density matrix [27, 28, 29]. 



2.6 Example in CFT: 2d massless Dirac fermion 

For our second example we consider the 2d massless Dirac fermion. This theory is sufficiently 
simple that, by borrowing heavily from existing results in the literature, we may compute 
the MI and TMI exactly. 



12 




° 1 



Figure 3: Depiction of the Euclidean path integral (A. 12). (I) shows the placement of 
delta function insertions used to construct the Euclidean analogue of Paxubx] (II) shows the 
placement for the Euclidean analogue of Paxub2- The twist fields (see §2.6) are also depicted, 
note that the orientation of the cuts at —(3/2 are reversed. 



Following the discussion of the previous section, our main task is to compute the Renyi 
moments tr^ c (pg.) starting from the expression for the Euclidean path integral for pc- We 
may use the replica trick to glue together path integrals for pc in order to obtain a single 
path integral for tr^ c (pg.) over an n-sheeted Riemann surface. For 2D CFTs we may use 
a further refinement of Calabrese and Cardy [4] who have shown that this path integral 
may be related to the correlation function of twist operators cr^(z) defined on the original 
space AIe- Let us denote by [u, v] a line on AIe with endpoints u, v and constant imaginary 
part, and let C be the union of N disjoint segments C = U^-Jitj, The trace tr% c (p£) is 
proportional to a 2iV-pt. correlation function of twist operators: 

ti nc (p c ) n oc (atMo5(vi) . . . a±(u N )a±(v N )) n . (2.41) 

The twist operators are conformal primaries with conformal weight 

12 V n J 12 n 

with c the central charge. The choice of which twist field is determined by the placement of 
the intervals. For the case considered by [4] where all intervals lie on the real axis (Imz = 0) 
the twist correlator is 

tr Hc {p c ) n oc (<t+(ui)<7~(ui) . . .a+(u N )cr~(v N )) n , Imia = Imvi = V % = 1, . . . , N. 

(2.43) 

For segments lying on Imz = —{3/2 the orientation of the twist fields is reversed. This is a 
consequence of the it prescriptions described in §2.5. 

In order to compute the MI and TMI we consider the three segments 

(2.44) 
(2.45) 
(2.46) 



A x = 


Wax- 


vai], 


Im u^i 


= ImvAx 


= 0, 


Bx = 


[ubx 


vbi], 


Ivcxubx 


= Im v B x 


= 0, 


B 2 = 


[U B 2 




lxxxu B 2 


= lmv B2 


= -P/2, 
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The Renyi moments for the associated density matrices are proportional to 

tr (PAi) °c (cr+( y u A i)(r^(vAi)) n , 

tr (Pfi 2 ) oc (a-~(ttB 2 )^(vfl2)> n . 



(2.47) 
(2.48) 
(2.49) 



These 2-pt. functions are completely determined by conformal symmetry and depend on the 
specific model only in the value of c which determines A n : 



Tie 



(3 sinh 

The Renyi moments of unions of the regions are proportional to 



ttQi-di) 

p 



2A„ 



HPaiubi) oc (<T+(uAi)<T n (v A1 )<T+(uBi)<r n (v B i))i 
tT (PAlUB2) oc <(T+(m A i) ( t-(w a1 )(T^(m B2 ) ( t+(i; B 2)>j 



(2.50) 



(2.51) 
(2.52) 



Since the 4-pt. functions are not fully determined by conformal symmetry the right-hand 
side of these proportionalities are model dependant. For the massless Dirac fermion 3 c = 1 
and the twist correlators are known [31] : 4 



+ {u l )a n {yi)a+{u2)cr n {v2)) { 



7r 2 e 2 sinh 


7r(«l— U2) 
P 


sinh 


7r(«l— V2) 
P 




f3 2 sinh 


P 


sinh 


tt(u\—V2) 

p 


sinh 


7r(ui— 112) 
P 


sinh 


TV(U2— V2) 

P 



2A r , 



, (2.54) 



[a+(ux)a n (vx)a n (u 2 )a+(v 2 )). 



7r 2 e 2 sinh 


7r(«l— U2) 


sinh 


7r(Dl— U2) 

P 




/3 2 sinh 


7r(ui— Ul) 



sinh 


7T(«1— U2) 

/9 


sinh 


■k(v\—V2) 

p 


sinh 


it(u2 — V2) 

p 



2A„ 



. (2.55) 



We can now quickly assemble these ingredients into the MI and TMI. The mutual infor- 
mation is given by 



MI(A 1 : B x ) = lim 



1 



1 1 - n 



log 



tr(p 



AlUBl, 



(2.56) 



3 As has been clarified in [30], we will be studying the unprojected Dirac fermion theory, not the modular 
invariant theory where fermion number is gauged. 



For this model the correlation functions of twist fields on the complex plane are given by 

,2A„ . - e 



(a+(ui)(r n (vi) . . . cr+ {u N )a n (v N )) p i an e = |detM|' 



My = 



(2.53) 



The correlation functions on the cylinder may be obtained by the conformal transformation w p i an e 
exp[2irZcyi/0\. 
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We know the Reyni moments up to constants of proportionality which are independent of (5. 
These constants provide an overall constant to the mutual information. The correct value 
of this constant is fixed by the requirement that the mutual information vanish as (3 — > 0. 
Inserting the expressions for the twist correlators we obtain 



MI(Ax : B x ) 



1 



log 



sinh 



tt(u a1 -u b1 ) 



sinh 



P 



sinh 



n(u A1 -v B1 ) 

P 



sinh 



n(v A1 -u B1 ) 

P 



(2.57) 



To clean this expression up consider without loss of generality the configuration uai < vai < 
ubi < Vbx'i we ma y then adopt the variables 



\uai - VAl\ 
and rewrite (2.57) as 



La, 



ubi - Vbi 



ubi - v A i = S, 



(2.58) 



MI(Ai:Bi) = -log 



For the TMI we similarly compute 

1 



sinh ^^a+s) ^ s j nn ^zl 



(L B +S) 
P 



sinh 



tt(L a +L b +S 



Mill) ( ?f 



(2.59) 



TMl(A 1 :B 2 ) = - log 



log 



sinh ( " ( " A1 ~" B2) 



cosh 



n(L A +L B +S) 



cosh ( ^ 



cosh cosh (*&j±&) 



(2.60) 



(2-61) 



In the second equality we have once again inserted the values (2.58), along with the fact 
that lmuB2 = Ito.vb2 = —0/2- Once again the overall constant is established by demanding 
that the TMI vanish as (3 — > 0. 

We provide representative plots of the MI and TMI in Fig. 4. As expected, the MI bounds 
the TMI from above. Both MI and TMI vanish as /3 — > 0; as (3 — > oo the MI approaches is 
maximum value while the TMI once again vanishes. 



3 Holographic t her mo- mutual information 

In this section we change gears and utilize the holographic techniques of the AdS/CFT 
correspondence [32] to study TMI from a different perspective. In particular, we use the 
holographic formula of Ryu and Takayanagi [6] to examine TMI for the 2d CFT dual to the 
non-rotating BTZ black hole in 3 bulk dimensions. 

The original RT formula [6] provides a holographic prescription for calculating entan- 
glement entropy in a strongly-coupled <i-dimensional CFT on a static spacetime Ai which 
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5 10 15 20 0.5 1.0 1.5 2.0 

P S 



Figure 4: Examples of the MI and TMI of the 2D massless Dirac fermion. Left: the MI (solid 
blue) and TMI (solid red) as a function of 0. In this plot La = Lb = 1 and S = 1/2. The 
solid black line denotes the maximum value of the MI. Right: the MI and TMI as function 
of separation S. Here La = Lb = 1 and f3 = 5. 

admits a dual description, in the usual sense of AdS/CFT, as a static d + 1-dimensional 
asymptotically AdS spacetime B = M. X M. Because both bulk and boundary admit a 
timelike Killing vector d t , we may use the canonical foliation into equal t surfaces T, t C B 
and <9Xi C Ai. Consider a region A C <9E t on the boundary. The RT proposal is that the 
entanglement entropy Sa is given by 5 

4(jr jv Ma 

where Ma C is a surface in the bulk geometry homologous to A. While this is only a 
conjectured formula, it has passed many nontrivial checks (see e.g. [33, 9, 34] for a review). 
When the system in question is no longer static (3.1) must be generalized to the covariant 
holographic entanglement formula of [11]. The clearest presentation of this generalization 
replaces minimal surfaces with saddlepoints of the area action for spacelike co-dimension 
two bulk surfaces homologous to A. In 3 bulk geometries these surfaces are simply spacelike 
geodesies, and the entanglement entropy is then again given by (3.1). 

In AdS/CFT thermal systems at high enough temperature 6 and without chemical po- 
tential are described by non-rotating eternal black holes. These manifolds have two exterior 

5 We are working in Einstein frame with the action normalized as S — 16 ^ Gn J y/—gR + • • ■ and ignoring 
higher-derivative corrections. 

6 It is well-known that AdS black holes at sufficiently low temperature ((3 > 2tt in 3 dimensions) have 
higher free energy than global AdS and therefore are thermodynamically unstable to a first-order phase 
transition [35] . We will restrict attention to systems at temperatures well above the transition temperature, 
i.e., /3 < 2?r. 
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Figure 5: The Penrose- Carter diagram for a maximally extended AdS-Schwarzschild black 
hole. The two dark shaded regions labeled R± are the regions covered by two AdS- 
Schwarzschild coordinate patches where the killing vector dt is timelike, bounded by the 
45° dashed lines indicating the horizon. The red jagged line is the curvature singularity, 
which bends inwards for AdS-Schwarzschild black holes in D > 3. 

regions outside the black hole horizon and correspondingly two conformal boundaries, see 
figure 5. While these manifolds do not admit a global timelike Killing vector field, they 
do admit Killing vector fields which are timelike in the exterior regions. The AdS/CFT 
interpretation in this setting was established by Maldacena [17] and in many way follows 
the classic work of Israel [36] . The two boundaries correspond to the two copies of the CFT 
in the TFD description of the thermal state. Due to the presence of timelike Killing vector 
fields on the boundaries the CFT is governed by a time-independent Hamiltonian. The en- 
tanglement entropy of the CFT state is realized geometrically as the entropy of the black 
hole. Basic properties of the TFD system described §2.2, such as the fact that all type-1 
and type-2 operators commute, are clearly realized through this bulk description. 

Several authors have used the RT formula to study mutual information in thermal systems 
holographically (see e.g., [11, 9]). In order to distinguish the holographic prediction from 
exact results in the field theory we refer to this as the holographic mutual information 
(HMI). Let dt be the timelike Killing vector field in a exterior region, and consider two 
regions A\,B\ C <9£t on the conformal boundary of this exterior region. The HMI between 
A\ and B\ is simply 

HMI(Ai : B,) := S A1 + S B1 - S A1UB1 , (3.2) 

with the holographic entanglement entropies computed from the RT formula (3.1). The 
entropies Sai and Sbi are given by VoI(Mai) and VoI(Mbi), respectively, where Mai,bi C E t 
are the minimal surfaces in the bulk homologous to A\ and B\. For the entropy Saiubi there 
are two candidate minimal surfaces: i) the disconnected surface Mai U M B i, and ii) the 
connected minimal surface M c homologous to A\ U B 2 . The entropy Saiubi is determined 
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by the lesser volume of these two choices. Thus we may compactly write the HMI as 

HMI(Ai : B x ) = — - max [Vol(M A1 ) + Vol(M B1 ) - Vol(M c ), 0] , (3.3) 

where the zero result is obtained when the Saiubi is given by the disconnected minimal 
surface. It is important to note that none of these surfaces cross the black hole horizon; the 
geometry behind the horizon is not needed for (or probed by) this computation. This is the 
analogue of the fact that mutual information may can be constructed from only one copy of 
the CFT - the second copy is traced out to reproduce the thermal density matrix. 

The holographic TMI may be computed in a similar manner to the HMI. Since thermo- 
mutual information involves operators of both type-1 and type-2, the holographic thermo- 
mutual information (HTMI) involves regions on each conformal boundary of the extended 
black hole spacetime. Extended black hole spacetimes do not admit a global timelike Killing 
vector field, so in general the entanglement entropies must be computed using the covariant 
form of the RT proposal. As already mentioned, however, in 3 bulk dimensions this does not 
affect the computation. Let Ai,B 2 G <9E t be two regions on different conformal boundaries. 
The HTMI between A\ and B 2 given by 

HTMI^x : B 2 ) = — — max [Vol(M A1 ) + Vol(M B1 ) - Vol(M c ), 0] , (3.4) 
4Gjv 

where now Mc is the the minimal connected surface connecting A\ to B 2 which in addition 
satisfies a regularity condition we describe momentarily. Because there is an Einstein- Rosen 
bridge connecting the two asymptotic regions, there do indeed exist co-dimension two sur- 
faces which connect A% to B 2 . In general, finding these minimal surfaces for maximally 
extended black holes analytically can be difficult; however, in 3 dimensions the computation 
is much simpler because BTZ black holes are quotients of AdS^. 

The regularity constraint we impose on Mc in the computation of HTMI is the Lorentzian 
analog of the homology constraint imposed in the RT formula. The homotopy constraint 
says that, for the geodesies corresponding to A\ U B 2 , there exists a spacelike co-dimension 
one surface r AluB2 , in the bulk such that dr = Ai U B 2 U Maiub2- This requirement that 
the surface is space like means that in AdS$ the ribbon between the geodesies can not 
twist around, as it would no longer be spacelike. This condition is crucial, as there are 
configurations where the geodesies corresponding to a twisted ribbon are shorter than the 
untwisted case (for instance, regions with angular width ~ ir on opposite sides of the circle) 
but are rejected due to failing the homology constraint. 

We are also not interested in surfaces which wrap the spatial S 1 multiple times - these 
are guaranteed to be longer than the untwisted "straight across" geodesic, and so can be 
ignored at leading order. 

3.1 BTZ geometry 

This section serves as a brief review of the BTZ geometry. We emphasize the description 
of BTZ in terms of an embedding space as this formulation makes our computations below 
quite simple. Similar treatments can be found in [37, 38]. 
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Recall that AdSs may be constructed as a 3-dimensional surface in a 4-dimensional 
embedding space IR 2,2 . Adopting the metric gAB = diag{— , +,+,—}, for this space, AdS^ is 
the universal cover of the hypersurface defined by 



-e = g AB x A x B = -(xy + (x 1 ) 2 + (x 2 ) 2 - (xy, (3.5) 

with £ the AdS radius. BTZ is a quotient of AdS 3 and also admits a simple description as an 
embedded surface. If we identify X 2 ± X 3 = e ±47r2///3 (X 2 ± X 3 ) in the embedding space R 2 ' 2 
then the BTZ geometry is given by the hypersurface (3.5). We obtain the induced metric 
on the BTZ geometry by introducing the dimensionless real coordinates (w+,w-,<f>) via 



X° 
X 1 



— w+ — W- 

1 + W+W- 

—w+ + w_ 

1 + W+W- 



1 + w+w. 

These coordinates have the respective ranges: 



X* = £l 1 - W+W - Isinh, 
1 + w+w. 

X 3 = £( ]~ w + w - ) cosh0 . (3. 6 ) 



4vr 2 

w+ 6l, W- E R, -1 < w+w- < 1, = + — . (3.7) 







The BTZ line element is 



ds 2 = -, rrr \-Adw+dw^ + (-1 + w + w^) 2 d(f) 2 ] . (3.8) 

(1 + w+w^Y 

Clearly w + and W- are null coordinates and is an angular coordinate. There is a conformal 
boundary at w + W- = —1, a horizon at w+W- = 0, and a conical singularity at w+w_ = +1. 
See Fig. 6. 

For our purposes it is more useful to employ AdS-Schwarzschild coordinates. Consider 
the regions R+ := {0 < w+, — 1 < w+W- < 0} and i?_ := < 0, — 1 < w + W- < 0} 

respectively; each region is outside the horizon and may be covered by a copy of the AdS- 
Schwarzschild chart 

1 + w + w- (3 

Z — Z HZ , Z fc [\J,Zh), Zh — —, 

1 — w+W- lix 

= z H (p, 9 = 9 + 27r, 

W + = ± V2 e'/-, = + f ^4^) V2 e-/-, (3.9) 

\z H + zJ \z H + zl 
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Figure 6: The Penrose-Carter diagram of BTZ. The darker shaded regions are R + and 
covered by the AdS-Schwarzschild patches (3.9). Note that unlike higher-dimensional black 
holes, the singularity is not bent in. The (■, •) indicate the signs of w + and W- in each region. 



1 



'-H 



dt 2 + 1 



-1 



dz 2 + de 2 



(3.10) 



The choice of sign in the last line of (3.9) is ± for R±. We can analytically continue from 
R + to R_ via t — > t — i/3/2. In these coordinates the conformal boundary is at z = and 
the horizon is at z = Zh- Using this normalization of time, we see that the horizon has 
a temperature The full BTZ geometry also has "upper" and "lower" wedges (recall 

Fig. 6); in these regions d t is spacelike. 

The extended manifold does not enjoy a global timelike isometry, but it does admit a 
boost isometry that is timelike in the regions R±; the boost vector is 



'1 + w + w. 
4w- 



"04 



(i + w + w^y 

4w+ 



zh 



(3.11) 



provides a 



which has opposite orientation in the two regions R±. Kruskal time T = W+J ^ w - 
global notion of time but 8t is not an isometry. 

As we mentioned previously, there is a first-order phase transition between BTZ and 
global AdS at a finite temperature. One can easily compare the free energies F = E — TS 
of the two solutions, 

n 2 £ 



AF 



BTZ 



AdS 



-f3- 2 + 



(3.12) 



2Gjv SGj\[ 

which means that BTZ is only the dominant saddlepoint for > -L. In what follows as 
we discuss (T)MI in holographic thermal systems we will therefore always be assuming we 
are in this high temperature regime. 

The regions near the conformal boundaries are most easily investigated in the AdS- 
Schwarzschild charts (3.10). Letting z = e and taking the limit e — > the other coordinates 
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behave as 



w + = ±{l-e/z H )e t ' ZH + 0(e 2 ), 
w_ = ±(-1 + e/z H )e-^ ZH + 0(e 2 ), 
l + w + w_ = 2e/z H + 0(e 2 ), 

1-w+w- = 2-2e/z H + 0(e 2 ), (3.13) 
where the ± sign is chosen in region R±. 

3.1.1 Embedding distance 

The embedding space endows the BTZ geometry with a convenient notion of distance. We 
define the embedding distance between two points on BTZ to be Q(X,Y) := gAsX A Y B j I 2 . 
The embedding distance is simply related to the length of the chord through the embedding 
space 

\\X -Y\\ 2 = -2£ 2 (Q(X,Y) + 1), (3.14) 
as well as the geodesic distance D(X, Y) on the BTZ manifold: 

D(X Y) — { ^ acosn t~ ®P^> ^)] spacelike separation , , 

' [ £acos [— Q(X, Y)] timelike separation. 

In the coordinates (3.8) Q(X, Y) may be written 

t^, v v , 2(w +l w_ 2 + w_ l w +2 ) + (1 - w +2 W- 2 ) cosh(0! - 02 ) , . 

9 Ai,A 2 = — r . (3.16) 

Both the embedding distance and geodesic distance generically diverge as one or both points 
approach the conformal boundaries. Using the limits (3.13) we compute the limit of the 
embedding distance when both points approach a boundary: 

Q(X U X 2 ) = QrCS(Xl ' X2) + finite, 
ei£2 



e rcg (x 1; x 2 ) :=4 



s(Xi,X 2 ) cosh ( — J — cosh 



zh J \ z h 



(3.17) 



where s(Xi,X 2 ) = +(— ) for points approaching the same (different) boundaries. For points 
which are spacelike separated in this limit the geodesic distance behaves as 

r 1 D(X 1 ,X 2 ) = - ln( ei e 2 ) + In [-26 rcg (X 1 , X 2 )\ + ..., X U X 2 spacelike, (3.18) 

with ellipses denoting terms that vanish as ei j2 — > 0. In examining these relations it is 
useful to recall that the AdS-Schwarzschild coordinate t flows in opposite directions on the 
two boundaries. For instance, under the action of a boost with positive rapidity both t± 
and t 2 increase, regardless of which region the points belong, while the difference t\ — t 2 is 
unaffected. 
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3.2 Holographic computation 



The embedding space description makes it quite easy to compute both the holographic 
mutual information and holographic thermo-mutual information. Both quantities may be 
computed from the dimensionless "cross ratio" 



J ■= r 1 [D(X Al ,X A2 ) + D(X Bl ,X B2 ) - D(X A1 ,X B1 ) - D(X A2 ,X B2 )] 



(3.19) 



which we define for four points X A y, X A2 , X B i, X B2 in R + Ui?_. To obtain the HMI let the 
four points be in R + with AdS-Schwarzschild coordinates X = (t, z, 9): 



X A \ = (0, e, 9 A i), X A2 = (0, e, 9 A2 ), X m = (0, e, 9m), X B2 = (0, e, 9 B2 ), 



and 



< 9 A1 < A2 < 9 B2 < 9 B1 < 2tt. 



The HMI is then given by 



HMI(Ai : B v 



AG 



max 



N 



Mm J, 

,<E->0 



(3.20) 
(3.21) 

(3.22) 



Taking the limit with the aid of (3.18) and (3.17) we obtain 

(l - cosh (l - cosh 



HMI(A! : By 



AG 



In max 



Mb 

zh 



cosh ^) 



cosh 



A8 2 \ 
zh I 



(3.23) 



where we have introduced the obvious notation 



A9 A = \9 A1 — 9 A2 \, A9 B = \9m — 9 B2 \, A9y = \9 A y — 9 B1 \, A9 2 = \9 A2 — 9 B2 \. (3.24) 

The ordering (3.21) is a convenient choice of labelling so that (3.23) is constructed from the 
correct choice of geodesies such that the surface Mc is homologous to A\ U B\. 

The HTMI is similarly obtained from (3.19) by considering, e.g., X A1 ,X A2 e R + and 
X B y,X B2 G R-. In their respective AdS-Schwarzschild charts the points may once again be 
labelled by (3.20) with angles satisfying (3.21). This choice of labelling selects the correct 
surface Mc for defining S A i uB2 - recall discussion in §3. The HTMI is then 



HTMI(Ai : B 2 ) :-- 
and taking the limit we obtain 

HTMI(Ai : B 2 ) = ^^lnmax 



AG 



max 



N 



\imJ(A : B), 

e-S-0 . 



(3.25) 



AG 



N 



cosh 

zh J 



cosh^ 

ZH 



1 + cosh ^ 

ZH 



1 + cosh ^ 

zh 



(3.26) 



As with the thermo-mutual information computed in the quantum mechanics and CFT 
examples above, the HTMI captures the basic properties of thermo-mutual information 
outlined in §2.3, though the agreement is not perfect. In particular we note: 



22 



14 



12 
10 

8 

6 

4 

2 



A</> 
S 

Figure 7: An example of the holographic mutual information (HMI) and thermo-mutual 
information (HTMI). For this plot the angular widths are A4>a = = A<p; S denotes the 
angular separation between the centers of A and B respectively. In blue are profiles of the 
HMI for various temperatures, with increasing temperature yielding a non-zero HMI at larger 
S. In red are profiles of the HTMI for the same temperatures, with increasing temperature 
yielding greater HTMI. Note that the HTMI vanishes when the 9 profiles of A and B have no 
overlap. For sufficiently low temperature, though still above the Hawking-Page transition, 
the HTMI vanishes for all S. These features are generic. 

i) By construction HTMI is non-negative. 

ii) From (3.26) it is clear that HTMI is bounded above for all configurations. The fact that 
HTMI is constructed from two regions on different boundaries makes it manifest that 
no UV divergences are encountered when A\ and B 2 are taken to have overlapping 9 
profiles. 

iii) Given two regions B\^ with equivalent 9 profiles, one region on each boundary, we 
obtain from (3.26) and (3.23) that HMI(Ai : B x ) - HTMI (At : B 2 ) > 0. In this sense 
HTMI bounds HMI from below. 

iv) We may readily verify from (3.26) that HTMI^! : B 2 ) tends to zero in the limit 
/3 — > oo. In fact, the HTMI vanishes for /3 ~ ^V9~a9b, which can occur well above the 
Hawking-Page temperature /3 = 2n. 

This last point, the fact that the HTMI can vanish at temperatures well above the Hawking- 
Page temperature, deserves further discussion. In §2.3 we made a point to note that at any 
finite temperature the TMI cannot vanish for all configurations. This follows from the fact 
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that the TMI bounds the correlations of type-1 and type-2 fields, and these correlations 
cannot everywhere vanish if Q is indeed entangled. The fact that the HTMI can vanish for 
typical regions at finite temperature is a large N artifact. Recall that the RT formula is 
believed to capture only the 0(N 2 ) contributions to the entanglement entropies. Therefore 
the HTMI is itself contains only the leading large N behavior of the TMI. Correlations 
between type-1 and type-2 fields are naturally O(N ) are therefore are technically zero at 
this order. We provide a representative plot of HTMI and HMI in Fig. 7. 

4 Discussion 

In this paper we have begun a study of holographic mutual information between regions on 
disconnected AdS boundaries. In AdS/CFT these bulk geometries are commonly interpreted 
as describing entangled states in tensor product (i.e., non- interacting) CFTs. Specifically, 
we have focused on spacetimes with two disconnected boundaries which describe entangled 
states in the thermo-field double formalism. The mutual information we compute holo- 
graphically is the thermo-mutual information (TMI) between 'physical' and thermo-double 
degrees of freedom. To better understand our holographic calculations we have examined 
TMI directly within field theory and have detailed a number of it's attributes in §2.3. We 
have also computed TMI explicitly in two examples: a quantum mechanical spin chain, and 
a 2d Dirac fermion. The holographic system is in remarkable agreement with these examples 
and the attributes listed in §2.3. Thus our study provides yet another positive test for the 
conjectured Ryu-Takayanagi formula [6]. 

Perhaps the biggest discrepancy between the holographic and exact TMI is the fact that 
the holographic TMI can vanish for any choice of region at finite temperature, contrary to 
our general expectation that TMI cannot in general vanish at finite temperature. This is due 
to the fact that the holographic computation captures only the leading behavior in N 2 , and 
so like the holographic mutual information exhibits features similar to a first order phase 
transition. Indeed, that it is precisely the same phenomena as for the holographic mutual 
information may be seen as follows. One may adopt Rindler coordinates on the boundary 
of global AdS and in doing so describe the zero-temperature pure state (from the global 
perspective) as a thermal state at the Rindler temperature. The CFT in one Rindler wedge 
may be interpreted in the TFD language as the thermo-double of the the CFT in the other 
wedge. In this set-up the thermo-mutual information is precisely the mutual information of 
the global perspective. 

As this is a very preliminary study of this type of entanglement in holographic systems, 
we close with a number of discussion points and future directions. 

4.1 Holographic TMI in higher dimensions 

A natural extension of this work would be to examine holographic TMI computed in higher 
dimensional black hole backgrounds. While conceptually a straight-forward generalization, 
the necessary computations are considerably more involved in d + 1 bulk dimensions with 



24 



d > 2. Partly this is due to the fact that in higher dimensions the covariant form of the RT 
formula is necessary and the extremal surfaces involved have dimension greater than two. 
In addition, in 2 + 1 bulk dimensions our computations were rather slick because of the fact 
that BTZ is a quotient of AdS-^. Regardless of these computational challenges, we expect 
similar physics in higher dimensions. Namely, we expect the holographically computed TMI 
to be non- vanishing even at leading order in N 2 , and we expect it to exhibit a first order 
phase transition at /3 ~ ttL where L is the length scale of the regions under consideration. 
The features of TMI listed in §2.3 are of course independent of dimension. 

4.2 Other asymptotically AdS-Schwarzschild geometries 

It has been pointed out by several authors [39, 40, 41, 42, 43] that in 2 + 1 Einstein gravity the 
maximally extended BTZ geometry is not the only solution which has one asymptotic region 
corresponding to that of the AdS-Schwarzschild chart (3.10). There also exist solutions 
with more than two asymptotic regions as well as those with only one asymptotic region 
but non-trivial topology behind the horizon. Geometries with more than two asymptotic 
regions may be thought of as describing different choices of purification for a system initially 
described by a thermal density matrix. These purifications enlarge the total Hilbert space 
to a space other than the TFD Hilbert space "Hi <g) H,2- We expect that basic features of 
these purifications may be determined by examining various HTMI-like objects. 

Solutions with only one asymptotic region but non-trivial topology are even more fasci- 
nating. The AdS / CFT intuition states that the CFT dual of these geometries must be in a 
pure state. However, the RT formula applied in any black hole spacetime predicts a mixed 
state: due to the presence of the black hole in the bulk, the RT formula computes for any 
region A that Sa ^ S^, indicating that the state is mixed. For these spacetimes it appears 
that the RT formula gives the incorrect result at leading order in N 2 , but this warrants a 
more careful analysis. 

4.3 Renormalized entanglement entropy and mutual information 

We have emphasized throughout our analysis that one important reason mutual information 
(and, for intersecting regions, the F-function) is an object of interest is that it is a manifestly 
UV finite quantity. Recently [44] has proposed a UV-finite "renormalized entanglement en- 
tropy" (REE). It is interesting to ask if the same physics contained the MI may also be 
recovered from linear combinations of the REE. For a 2d system [44] define the REE for a 
subregion A with length scale L to be Sa = L8lSa- The naive "renormalized mutual infor- 
mation" is then L8l{Sa + Sb — Saub)- However, this quantity does not appear to capture 
the same physics; in particular, the derivative with respect to L annihilates the dependence 
on the scale-invariant cross-ratio. For zero temperature CFTs the mutual information is a 
function only of this cross-ratio, so this renormalized mutual information vanishes. For CFTs 
at finite temperature the naive renormalized mutual information is non-zero but again does 
not depend on the cross-ratio and so does not contain information about correlations, only 
information about the geometries under consideration. We mention this here in the hope 



25 



that it spurs a better understanding of the relationship between the proposed REE and the 
UV-finite quantities considered in our analysis. 
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A Euclidean path integrals for reduced density matri- 
ces in TFD 

In this appendix we derive the results quoted in §2.5. We begin by reviewing the relationship 
between the TFD path integral and the Euclidean path integral. Boundary conditions must 
be supplied in order for any path integral expression to be well-defined; for thermal systems 
the relevant boundary condition is the KMS condition [45] . Denoting the time translation of 
4> by 4>(t) = e lHt (p(0)e~ lHt and suppressing other spacetime arguments, one readily deduces 
from (2.10) that correlators of the thermal state Q at inverse temperature satisfy 

(<Kt)jd a =<x<Kt + iP))n, (A.l) 

for two operators (f>(x) and x( x )- A consequence of (A.l) is that correlators with respect 
to Q are boundary values of analytic functions in the complex time domain Imt e (0, — 0) 
satisfying the condition (A.l). Therefore we may seek a path integral generating function 
Zc[J] for the path-ordered correlation functions (P(j)(xi) . . . <p(x n )) Q . The 'path' refers to the 
time integration which is traversed from Imt = to Imt = — (3 along a contour C; a more 
careful examination of the analyticity properties of such correlators reveals that the path 
must have monotonically non-increasing Imt [24]. The path integral may be constructed in 
the usual manner: 



Z C [J] := Z- L / [V<j>] exp 



c 



iS c [<j>]+i dt d D - 1 x^-g{x)J{x)4>(x) 



(A.2) 



D is the number of spacetime dimensions and Z = Z[0]. The subscript Sc[</>] denotes that, 
like the source term written explicitly in (A.2), the time integration in the action is over 
the path C. The boundary conditions imposed on the field integration is 0(t) = 0(t + i/3). 
Path-ordered correlation functions may be computed by taking functional derivatives with 
respect to the source J(x): 



iSJ(xi) ■ ■ ■ i5J{x n ) 



(A.3) 

j=o 



The two choices of time contour relevant for our discussion are shown in Fig. 2. The 
first choice is the Euclidean contour which is traversed along the imaginary t axis — > 
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This path integral determines the Euclidean correlators (<f>(xi) . . . <f)(x n )) E which define the 
Euclidean 'state.' Explicitly, 



Z E [J] := Z- 1 j [V4>] exp -S E [4>] - j dr j d D - x xyfg(£)J(x)<j){x) 



(A.4) 



A second time contour is the TFD contour - see Fig. 2. With ti and tf real initial and final 
times respectively, this contour is traversed: i) U — > tf, ii) tf tf — i(3/2, iii) tf — i/3/2 — > 
t,i — i/3/2, iv) ti — i(3/2 — y ti — i/3. Under reasonable circumstances (see discussion §2.4 of 
[24]) the path integral factorizes into two path integrals, one involving the horizontal legs of 
the contour, the other involving the vertical legs. If one is interested only in correlators of 
operators lying on the horizontal contours then the latter path integral is simply 1 due to 
the overall normalization. The path integral for the horizontal segments is 



ZtfdIJiiM ■= Z T p D J [Vfa^Vfa] exp 



+ iS[<p 1 }-iS[(t )2 ] 



+ i / d D x^/g{^{Jl{x)(j) 1 {x) - J 2 {x)(j)2{x)) 



(A.5) 



Here (j)\{x) denote fields on the real axis while 4> 2 {x) fields have Imt = —/3/2. As our 
notation implies, these fields correspond to the type-1 and type-2 operators described in 
§2.2. Correlators of the operators are generated via: 



[[T<j>i(xi) . . . <f>i(x n )] [T(f) 2 (yi) . . . MVm)] > n 

5 n+m Z TFD [J u J 2 ] 



iSJ^x) ■ ■ ■ i5J 1 {x n )(-i)5J 2 {y l ) ■ ■ ■ (-i)5J 2 (y m ) 



(A.6) 



Ji=J 2 =0 



We may now construct path integral representations of a reduced density matrix in the 
TFD system. Let us first consider the density operator p^i = p(Q', Ai(A)). Without loss 
of generality we let the region A C So- For a basis in Wax we use the set of states {ipAi} 
specified by their field configuration of type-1 fields ip\(x) on A: 



: A1 



{4>\(x) = if)i(x) I x e A}. 



(A.7) 



We denote the matrix elements of pai in this basis PaiO^'aiU'ai)- Following [3], the path in- 
tegral representation of Pai^'aiI^ai) is then easily constructed by inserting delta functionals 
into the path integral (A.5): 



PAi{ip A i\i>Ai) 



r-l 



[V<j> x ][V<t> 2 ]{ expl+iSfal-iSlfa]] 



J! 5(0i (0 + ) - ^(O + ))5(0 1 (O-) - ^i(o-)) 



.xeA 



(A. 
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We have suppressed the dependence of spatial coordinates. The limits 0^ denote how t 
should be taken to approach t — 0, and reflect the fact that we define ip' (ip) in the limit 
t + ^- ) . The analogous expression for the density matrix p^ 2 , in an analogous basis {4>A2}, 
is 

PA2^'a2\M = Z- 1 J [VfaWVM exp - iS[(j> 2 }} 



J] 6(<h(-iW + + ) - $(-tf/2 + 0+)) 
<f(&(-i/9/2 + 0") - V2H/V2 + 0-)) 



(A.9) 



The density matrices paiubi and paii>B2 are constructed from the obvious generalizations of 
these expressions. 

We would like to compute these path integrals in Euclidean signature; this requires that 
we deform the TFD time contour to the Euclidean contour. To do so we recast the t + 0^ 
limits in the expressions above in terms of path ordering rather than real time ordering. Let 
X(t) denote a function that decreases monotonically along the path from ti to U — i(3. By 
rewriting the TFD path integral in terms of X(t) it is then straightforward to deform the 
time contour to the Euclidean contour. We then let A = Imt. For instance, after performing 
this step we obtain a Euclidean expression for the density matrix p A1 : 



= Z' 1 J[D<f>] 



-S E [ 



n 5(0(0-) -^(o-))5(0(o + )-^(o + )) 



(A.10) 



where now we label only the dependence on Imt. By the same steps we obtain the analogue 

Of PA2- 



JazWaWa) = Z- 1 J[V4>]e~ s ^ 



n sM-p/2 + o + ) - v 2 (-/?/2 + o + )) 

x£A 

5(0(-/3/2 + O-)-^ 2 (-/3/2 + O-)) 



(A.ll) 



Notice that the ie prescriptions defining the states have switched. Finally, we record the 
formula for the Euclidean analogue of the density matrix Paiub2- This is the key formula 
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necessary to compute the thermo-mutual information via the replica trick: 
Z- ] I [£>0]e- 5E[ ' 



nwoi-^(o-)w(o + )-# + )) 

xeA 

H J(0(-/5/2 + 0+) - ^(-/3/2 + 0+)) 



5(0(-/3/2 + O-)-V(-/3/2 + O-)) 



(A.12) 
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